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We describe the dissociation of a diatomic Feshbach molecule due to a time-varying external 
magnetic field in a realistic trap and guide setting. An analytic expression for the asymptotic state of 
the two ultracold atoms is derived, which can serve as a basis for the analysis of dissociation protocols 
to generate motionally entangled states. For instance, the gradual dissociation by sequences of 
magnetic field pulses may delocalize the atoms into macroscopically distinct wave packets, whose 
motional entanglement can be addressed interferometrically. The established relation between the 
applied magnetic field pulse and the generated dissociation state reveals that square-shaped magnetic 
field pulses minimize the momentum spread of the atoms. This is required to control the detrimental 
influence of dispersion in a recently proposed experiment to perform a Bell test in the motion of the 
two atoms [C. Gneiting and K. Hornberger, Phys. Rev. Lett. 101, 260503 (2008)]. 
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I. INTRODUCTION 

The emerging field of ultracold atoms makes it pos- 
sible to perform hitherto unprecedented experiments on 
the quantum nature of material objects, thus introduc- 
ing a new quality compared to previous quantum ex- 
periments with immaterial photons. This is owing to 
the experimental control over the atomic state reaching 
a level where quantum mechanical effects become rele- 
vant. Given that Bose-Einstein condensates (BECs) are 
nowadays routinely produced, they can serve as an ideal 
starting point, for example, to probe condensed matter 
physics in a highly controlled environment provided by 
optical lattices QJII]. A considerable extension of the 
scope of such experiments has been achieved by exploit- 
ing Feshbach resonances to produce molecules and molec- 
ular BECs (mBECs) [SHPIj. The ability to control the 
interaction between the atoms with an external mag- 
netic field permits one to realize, for example, the BCS- 
BEC crossover [1] or Efimov states [13], and to establish 
coherent atom- molecule oscillations [THUS] , hinting at 
quantum-coherent chemistry. 

Beyond that, it has been recognized that the con- 
trolled dissociation of mBECs can also serve as a re- 
source for entangled atom pairs that would permit the 
demonstration of nonclassical correlations on a macro- 
scopic level [EH . The use of such a controlled disso- 
ciation for spectroscopic purposes has already been dis- 
cussed theoretically [HU [52] and demonstrated experi- 
mentally .23, .24]. In a recent proposal, we investigated 
the possibility of using an arranged dissociation of Fes- 
hbach molecules in order to violate a Bell inequality in 
the motion of two atoms [25]. There, a sequence of two 
short magnetic field pulses dissociates a single molecule 
out of a dilute mBEC such that each atom is delocalized 
into two macroscopically distinct wave packets propagat- 
ing along the laser guide. The associated dissociation- 
time entangled (DTE) state can be considered the mat- 
ter wave analog of a Bell state whose capability to 



yield nonclassical correlations can be revealed by an in- 
terferometric protocol reuniting the wave packets on each 
side [57J 121] ■ The violation of a Bell inequality, however, 
imposes stringent requirements on the DTE state, which 
makes it necessary to know the generated dissociation 
state in detail [26] . 

The theory of Feshbach molecules is mainly concerned 
with an accurate description of the molecules in the 
bound regime and its vicinity, where the atoms interact 
strongly. For instance, it has been investigated in great 
detail that the molecular state at constant magnetic field 
exhibits universal properties in the vicinity of a Fesh- 
bach resonance, as well as the association and dissocia- 
tion behaviors under a linear magnetic field sweep, with 
emphasis on the converted fraction (see [HI 130] and 
Refs. therein). An appropriate description of the situa- 
tion proposed in [25] , however, requires a detailed knowl- 
edge of the dissociation state for a more general tempo- 
ral behavior of the magnetic field, varying on short time 
scales compared to the inverse resonance width. While 
in the interaction regime of the atoms one is then forced 
to resort to numerical methods [31] , the restriction to the 
asymptotic situation (i.e. for large interatomic distances 
and times long after the dissociation process) comes with 
significant simplifications that permit one to evaluate the 
generated dissociation state analytically. 

In this article we present a coupled-channels formu- 
lation of the dissociation process of an initially trapped 
Feshbach molecule exposed to a time-varying homoge- 
neous magnetic field. These single-molecule dynamics 
yield an adequate description in the relevant case of the 
dissociation of a single molecule out of a dilute mBEC, 
where interactions between the molecules and statisti- 
cal effects do not play a role. The formulation includes 
the center of mass motion required for the complete de- 
scription of the two-particle system. The dissociation is 
considered to take place in a shallow trap from where 
the dissociated atoms are injected into a strong guiding 
potential, which confines them to the longitudinal mo- 
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tion along the guide axis. Based on the time-dependent 
coupled-channel equations, we derive an analytic expres- 
sion for the asymptotic dissociation state of the ultracold 
atoms, allowing us to investigate the relation between the 
form of the applied magnetic field pulse and the resulting 
two-particle dissociation state. While Gaussian-shaped 
dissociation pulses turn out to yield a rather bulky mo- 
mentum distribution of the two-particle state, we demon- 
strate that idealized square-shaped magnetic field pulses 
optimize the momentum distribution with respect to its 
sharpness. 

The structure of this article is as follows. In Section 
[IT] we describe the assumed geometry of trap and guiding 
potentials, which is the natural configuration for dissoci- 
ation experiments. The time-dependent coupled-channel 
equations are formulated in Section[TTTJ and then reduced 
to an integro-differential equation for the closed-channel 
amplitude and an associated equation for the background 
channel state. Section [TV] shows how the asymptotic dis- 
sociation state can be expressed in terms of the Fourier 
transform of the closed-channel amplitude. The approx- 
imate dynamics of the latter for a given shape of the 
magnetic field pulse is described in Section |Vj this al- 
lows us, in Section [Vl] to discuss the optimal form of the 
pulse if a narrow momentum distribution of the disso- 
ciation products is required, such as for the test of the 
Bell inequality proposed in [35] . Applications of the Fesh- 
bach dissociation scheme as a resource of entangled atom 
states beyond the scenario in [35] are discussed in the 
Conclusions. 



II. TRAP AND GUIDE CONFIGURATION 

Before presenting the coupled-channels formalism, we 
outline an experimental setup that meets the conditions 
assumed in our dissociation scenario. An explicit quan- 
titative elaboration was given in |32) . 

Consider a BEC of Feshbach molecules in an optical 
dipole trap, which can be established by two perpendic- 
ularly crossing laser beams, see Figure [I] a). The weak 
trap laser guarantees longitudinal confinement within the 
wave guide produced by the strong guiding laser. We take 
the BEC, of the order of 10 2 molecules, to be sufficiently 
dilute such that one may neglect interactions between 
different molecules. This may be accomplished by choos- 
ing a mBEC of fermionic constituents, for example 6 Li2, 
where Pauli blocking further reduces the effect of inter- 
molecular interactions and thus enhances the lifetime of 
the mBEC [S3]. 

The molecules can then be considered to be in a prod- 
uct state with the center of mass motion given by the 
ground state of the trap, whereas the relative motion is 
in a bound molecular state. The latter can be turned into 
a Feshbach resonance by varying the external magnetic 
field, allowing one to dissociate the atoms in a controlled 
way. By applying one or several appropriately chosen 
dissociation pulses, a single molecule dissociates into two 
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FIG. 1: (color online) a) Setup for generating pairs of mo- 
tionally entangled atoms by the Feshbach dissociation of 
molecules. Initially, a BEC of on the order of 10 2 Feshbach 
molecules resides in a dipole trap produced by two crossing 
laser beams. An externally controlled, homogeneous mag- 
netic field induces the dissociation of one molecule per trial 
and thus generates a pair of atoms moving along the laser 
guide at a velocity on the order of 1 cm/s. The asymptotic 
two-atom state in the wave guide is determined by the trap 
and guide geometry and by the shape of the dissociation pulse, 
b) The dissociation pulse promotes the trapped molecule to a 
pair of counter-propagating atoms. The laser parameters are 
chosen such that the energy supply from the magnetic field 
sweep exceeds the trap laser potential, while it cannot lead 
to an additional transversal excitation if the kinetic energy 
of the atoms _Bkin remains below the harmonic energy gap, 

fttJG > Skin- 

counter-propagating atoms on average, and we postselect 
the single-dissociation events. 

The pulses provide the atoms with a kinetic energy 
sufficiently large to overcome the trap potential in longi- 
tudinal direction, but still below the threshold to get be- 
yond the ground state transversally, see Figure [T]b). This 
way we may end up with two dissociated atoms, counter- 
propagating with a velocity on the order of 1 cm/s along 
the guiding laser axis, whose two-particle state is deter- 
mined by the initial state of the molecule and the dissoci- 
ation pulse shape. By applying sequences of dissociation 
pulses, one can design highly nonclassical, motionally en- 
tangled states that can further be processed for various 
fundamental tests of quantum mechanics, such as to vi- 
olate a Bell inequality [25 . 



III. COUPLED-CHANNELS FORMULATION 

The dynamics of Feshbach molecules is appropriately 
described by the coupled channels formulation. In our 
case, the channels are characterized by different nuclear 
spin configurations of the two-atom system. The different 
magnetic field dependences of the corresponding energy 
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levels make it possible to manipulate the system exter- 
nally. In the following we adopt the notation and the 
conventions from [33] as far as possible. As a novelty, 
we must include also the center of mass motion of the 
two atoms. We assume that the magnetic field remains 
always in the vicinity of a single Feshbach resonance, al- 
lowing us to restrict the description to two channels: the 
closed channel of the energetically more favorable spin 
configuration supporting the molecular bound state, and 
the background channel, where the dissociated atoms are 
asymptotically free. 



A. Hamiltonian and coupled-channel equations 

The two-channel Hamiltonian for the two atoms (here 
taken to be of equal mass m) can be written as 

H tot = ff c i|cl)(cl|+if bg |bg)(bg| 

+W|bg)(cl| + ^|cl)(bg|, (1) 

with the closed-channel Hamiltonian 

Ha = -^Vi-^V 2 2 +V c} (\xi-x 2 \,B(t)) 

+V T (x 1 ) + V T (x 2 ) + V G (x 1 ) + V G (x 2 ) (2) 

and the background (open) channel Hamiltonian 

& 2 & 2 

Hhs = ~2^ V '~ ^ V 2 + V bg (\x 1 -x 2 \) 

+V T (x 1 ) + V T (x 2 ) + Vb(xi) + V G (x 2 ). (3) 

Here Vr(xi) and V G {xi) denote the trapping and guid- 
ing laser potential, respectively (V G may contain a lin- 
ear shift due to the gravitational potential). Vd(|a3i — 
x 2 \,B(t)) and Vb g (|a;i — x 2 \), on the other hand, describe 
the interatomic potentials for each channel. In general, 
these potentials differ for the different channels, reflecting 
their dependence on the spin configuration. In the chosen 
convention, the zero of total energy is defined in absence 
of the laser potentials by the background channel dissoci- 
ation threshold with the center of mass at rest. Then only 
the closed-channel potential V^idaJi — x 2 \, B(t)) depends 
on the external magnetic field B(t), describing an overall 
shift with respect to the background channel dissociation 
threshold. 

The off-diagonal elements W denote the energies asso- 
ciated with the spin exchange interaction and provide the 
inter-channel coupling. We assume W to be diagonal in 
position (i.e. independent of momentum) and to depend 
only on the interatomic distance \x\ — x 2 \ from now on. 
For the following it is useful to reformulate the Hamilto- 
nian in center of mass (cm) and relative (rel) coordinates, 
x cm = (xi +x 2 )/2 and a; ro i = Xi — x 2 , respectively, with 
total mass M — 2m and reduced mass fi = m/2. The 



closed-channel Hamiltonian thus reads as 

H c i = -|^VL - ^V? cl + V cl (\xUB(t)) 

+V T (x cm + x 1 . c i/2) + V T ( /2) (4) 

+V G (x cni + x Te i/2) + V G {x cm - x Tel /2), 

and similar for H hg with V c i(\x re \\, B(t)) replaced by 
Hgd^rciD- Note that the center of mass motion is only 
indirectly affected by the homogeneous external magnetic 
field -B, due to the presence of the trapping potentials. 
Writing 

|*tot(*)) = ®cl(x c m,X re \,t)\ Cl) + $bg(x cnl , X Ieh t) | bg) , 

we can infer from the time-dependent Schrodinger equa- 
tion the coupled channel equations 

ift|$d(*)) = HjB(t))\**{t)) +W\$ hg (t)) (6) 
ifi|$ bg (i)) = H hg \$ hs (t)} + W\$ d (t)}. 

The closed channel and background channel two-particle 
state components |$ c i(t)) and |$bg(i)) are thus not nor- 
malized in general. It is reasonable to assume the trans- 
verse motion of the atoms to be constrained to the region 
of validity of the harmonic approximation of the guiding 
laser potential, and the harmonic approximation to be 
applicable also for the trap laser potential in the closed 
channel. We then can write (adopting cylindrical coor- 
dinates) 

V G { Xl ) + V G {x 2 ) = -U 0iG + — u G p\ - U 0iG + —uJ G p 2 2 

TTt TTL 

V T (xi) + V T (x 2 ) = -f/ ,T + y oJ^zl - f7 ,T + —u\z%. 

(7) 

The harmonic approximation comes with the virtue not 
to couple the center of mass and relative motion, 

Vfc(ssi) + Vb(a:a) = -2 C/ ,G + ^gpL + f ^p r 2 el 

V T { Xl ) + V T {x 2 ) = -2 U 0<T + z c 2 m + 14 4j. 

(8) 

The initial bound state can therefore be taken to be sep- 
arable with respect to its center of mass and relative mo- 
tion. Of course, the harmonic approximation for the trap 
laser potential breaks down when the dissociated atoms 
leave the trap. Then, the center of mass motion ceases 
to be bound by the trap potential but rather undergoes 
a free propagation resulting in a dispersive broadening 
on the spot. So, even though initially only the relative 
motion is affected by the external magnetic field, its ef- 
fective coupling in the background channel to the center 
of mass also couples the motion of the latter indirectly 
to the external magnetic field. 
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B. Single-resonance approximation 

It is legitimate [29] to take the relative motion of the 
closed channel state component to be proportional to the 
underlying bare resonance state |0res)j which is defined 
by 



2/i 



V r 2 cl + V c ,(|aw|,B(t)) 



0res (*^rcl) 

E ies (B(t))<t> iea (x iel ). (9) 



Note that the time-dependent external magnetic held af- 
fects only its energy, which is taken to vanish at the res- 
onance, E lea (B res ) = 0, such that E ICS describes the en- 
ergetic offset of |0 rcs ) from the background channel dis- 
sociation threshold. 

Since the laser potentials vary weakly over the spatial 
extent of |0 res ), this resonance state remains a valid ap- 
proximation even in the presence of the trap. We assume 
that |^> r es) is spherically symmetric and thus supports an 
s-wave resonance. If we further take into account that 
the center of mass motion of the closed channel is com- 
pletely determined by the longitudinal and transversal 
trap ground states \iPt) and l^o™), respectively, we can 
write 



where 



(10) 



fi2 n2 

V 

2M cm 



+ 2V T {z cm ) + 2V G (p cm )]V> T (2cmM^(Pcm) 

2U ,t + huj T /2 - 2Uo. G + fkJ G ]tp T (z cm )(poo{p cm ). 

(11) 



In the single-resonance approximation (extended by the 
trapped center of mass motion) |l0| ) , the spatial shape of 
the closed channel state component is not affected by the 
external magnetic held and therefore time independent. 
The closed-channel amplitude C(t) therefore captures the 
complete effect of the time- varying magnetic held. Using 
the single-resonance approximation ( |10[ ) and introducing 
the abbreviation U c \ — — 2[/o,t + h^Jr/2 — 2[/o,g + h^G, 
we can thus rewrite the coupled-channels equations Q 
as 



ihdtClt) = [E Ies {B(t)) + U d ]C(t) 

+(^m m m es \w\$ hs (t)) 



(12) 



(ihd t -H hs )\^ hs (t)) = C(t)W|^)bSo)|^™>. (13) 



function of the background channel Gh s (t, t'), which sat- 
ishes 



(ihd t -H hg )G hg {t,t') = 6(t-t'). 



(14) 



If we further make use of the connection between the 
(retarded) Green's function and the time evolution op- 
erator, Gbg(M') = Ub g (t, t')Q(t — t')/(iH), we can write 
the background channel state component as 



I'M*)) = 



i 



di/c(i/)u hs (t^)w\^T)m\<t>, 



(15) 

In our scenario, the boundary conditions prohibit a ho- 
mogeneous solution. Physically, this reflects the fact that 
the closed channel is the only source for the background 
channel, in particular there are no further sources at in- 
finity (for example incoming and scattered particles). A 
closed equation for the closed-channel amplitude C(t) 
arises from inserting the formal solution (15) into ( |12[ ), 
which yields 

[iW t - E Ies (B(t)) - U cl ]C(t) = ^ f dt'C(t') 

X (^ T |(^|(0rc S |W bg (t,i')^I^T)|^)|0rc S >. (16) 



With ( 15 1 and (|16[) we arrived at a decoupled set of equa- 



tions that divides the determination of the background 
channel dissociation state into two parts, first solving 



(16) for C(t), then using the solution in (15) for the cal- 
culation of |<&b g ) • The dynamics of the closed-channel 
amplitude C(t) is explicitly driven by the external mag- 
netic field B(t), reflected in the left-hand side of (16). 



The bare background channel Hamiltonian -Hbg, on the 
other hand, does not depend on the external magnetic 
held, which will allow us in Section |IV| to expand the 
right-hand side of ( 15 1 in terms of its (time- and coupling- 
independent) energy eigenfunctions. The right-hand side 
of (|T6j) describes the back action of the background chan- 
nel state component on the dynamics of C (t) due to the 
coupling W. A solution to ( 16 ) will be given in Section 

m 

As a last remark, we note that it might seem sugges- 
tive to first solve the time-independent coupled-channel 
equations for a stationary magnetic field B, and then to 
take the corresponding static decay rate to describe the 
dissociation in the time- dependent case B(t), as it was 
done in [25] . However, we will see later in this article 
that this quasi-stationary approach is not sufficient for 
our purposes. 



IV. ASYMPTOTIC DISSOCIATION STATE 



Formal Green's solution 



Interpreting the right-hand side of ( 13 1 as a source 
term for the background channel state component 



®bg{t)} suggests to solve (13) formally using the Green's 



The formal expression (15) describes the background 
channel state component in full generality. In the sce- 
nario described in Section [TTJ however, we only need 
to know the dissociation state for large interatomic dis- 
tances and for times long after the dissociation process. 
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Moreover, the relevant dissociation states are sharply 
peaked in the ultracold regime, in the sense that the 
width of the momentum distribution is much smaller 
than its average momentum, because only then they are 
useful with respect to further employment such as the 
Bell test in the motion [25 . The above restrictions ad- 
mit significant simplifications that permit us to provide 
an analytic expression for the dissociation state in the 
asymptotic regime. 



with the initial state in energy representation given by 
C(E/h)(E\ hg W\ip T )\^ofi)\<l>res)- We will see that this ex- 
pression describes two counter-propagating atoms with 
well-defined momenta if C(lu) is peaked at an energy in 
the ultracold regime. The corresponding wave functions 
then have de Broglie wavelengths and spatial extensions 
on the order of micro- to millimeters, all being features 
desired for applications such as further interferometric 
manipulation. 



A. Large time limit 



B. Quantum numbers 



As a first step, we expand the background channel 
time evolution operator U\, g {t,t') in an appropriate en- 
ergy eigenbasis of the background channel Hamiltonian. 



Note that the energy eigenvalues in (171 are highly 



U hg {t,t')=c 



-iffbg(t-t') _ ~-iE(t-t')/h 



= 2. e 

E 



l \E) hg (E\ hg , 



(17) 

where -ffb g |-E)b g = E\E) bg ; adequate quantum numbers 
for our setup will be specified later in this article. Note 
that the involved vectors are two-particle states. As men- 
tioned above, this representation is only possible for the 
bare background Hamiltonian, which does not depend on 
the external magnetic field. Since at large interatomic 
distances only the continuum states survive, we can drop 
the bound states in (15) for large |a? re i| and get 



(x cm ,x Tel \$ hs (t)) ~ I I dt ' C '( t ') ( 18 ) 

xe -iE(t-t')/h {Xc ^ Xi ^ E)hs 

x (£| bg ^ T ) raises). 

The sum over the energy eigenstates takes into account 
that the zero of energy is defined in absence of the confin- 
ing lasers, whose presence shifts the (longitudinal) con- 
tinuum threshold by an offset £7bg = — 2t/o,G + 2Hujq. 
We arrange the dissociation such that after its comple- 
tion the magnetic field persists at a base value Bq below 
the dissociation threshold. The closed-channel amplitude 
then shows a simple time dependence in the large time 
regime, C(t) = Coex-p(—iEot/h), such that for times long 
after the completion we can replace the upper integration 
boundary by infinity without modifying the integral for 
energies above the dissociation threshold. This allows us 
to interpret the integration over t' as the Fourier trans- 
form C(oj) = f_ dtexp(iu)t)C(t) of the closed-channel 
amplitude C{t), yielding 

(x cm ,x Iel \$ he (t)) « i V e- iEt ' h C{E/h) 

aSrel — >00 III * ' 

x (E\ hs W\^ T )\^)\ct> reB ) 

\E) hg . (19) 

We thus find that the asymptotic dissociation state can 
be interpreted as evolving in the background channel 



degenerate; in order to single out a unique energy ba- 
sis, we choose as asymptotically well-defined, commuting 
observables the complete set p cm , p rc \ and H±, where 
H±_ denotes the transversal Hamiltonian. Assuming har- 
monic transversal confinement, we can write 



E i s >b g (£ibg 



E>U he 



E 



G '"'G 



dp c 



dp 



i'd 



x \Pcim n Q ' m G i-Prcb^G ' m G )bg 
x (Pcm , n-Q 11 , m G m , p re i , rag 1 , mg 1 1 bg . 



(20) 



Here hq = 0, 1, . . . denotes the radial occupation num- 
ber and tuq = —hq, . . . , uq the corresponding angular 
momentum. Since we assume that W|<^ res ) is spheri- 
cally symmetric (s-wave resonance), W|^T)|^oo)l0res) 
is both cylindrically symmetric and symmetric under 
exchange of the two particles, hence only the states 
\p cm ,n,0,p le i,n,0) sharing these symmetries can be oc- 
cupied. 

If we take C{E/K) to be peaked at a sufficiently low 
energy and if the laser guide is chosen appropriately, we 
can ensure that only the transversal ground state is ener- 
getically accessible. In this single-mode regime we have 

{x cla ,x Iel \$ bg (t)) 

J — oo J — oo 



|x rol |->oo I 
t—t-oo 



x C 



Pc 2 m/2M+^ ol /2^ \ 

h ) 



x (p cm , 0, 0,p rcl , 0, 0\ hg W \ip T ) |<^™}|</>r 

X (^cm ) *^rel |Pcm,0,0,p re l,0,0)bg. 



(21) 



Even taking C(E/h) to be sharply peaked in energy, it 
still admits the whole class of degenerate eigenstates that 
fall into that energetic region. As we show next, the 
matrix element (p cm , , , p ro i , , 1 bg W \ ip T ) | V^o™ ) I <^rc S ) ef- 
fects a further restriction of the accessible eigenstates, 
yielding the physically appropriate description of the sit- 
uation. 
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C. Asymptotic center of mass motion 

The main effect of the iongitudinal trap is to correlate 
the motion of the interatomic distance with the center of 
mass, since the center of mass evolution is confined for 
small interatomic distances, while it is free for distances 
beyond the size of the trap. For a sufficiently narrow 
energy distribution this merely results in a global time 
delay for the start of the free center of mass propaga- 
tion, which may be neglected at the time scales of the 
asymptotic regime. The longitudinal center of mass mo- 
tion can thus be described by the momentum eigenstates 
|pcm), that is by 



|p cm ,0,0,p rc i,0,0) bg « |Pcm}|^o)brel,0,0) 



bg- 



Here |p re i, 0, 0)b g denotes the eigenstates of the Hamilto- 
nian for the relative motion, which still contains the traps 
and the interatomic potential Vb g . The dissociation state 
then reads 



(x cm ,x Ie i\$ bs (t)) 



\x te \\—>oo ih 
t— J-oo 



OO /"OG 

dpcm / dp rol 



x e -i(^bg+pL/2M+p r 2 ol /2p)t/S 
X (Pcm|V , T)(Prel 1 O,O| b gW'|0 reB ) 

x (z cm 

|Pcm)(Pcm|yo!o> 

^ {^rel 

|p re i,0,0) bg . (22) 

The matrix element (pcmlV'T), given by the longitudinal 
harmonic trap ground state, guarantees that the center of 
mass motion remains centered at vanishing momentum, 
as required by momentum conservation in the dissocia- 
tion process. Since we take the energy to be in the ul- 
tracold regime of the background channel potential Vb g , 
the matrix element (p re i, 0, 0|b g VF|</>res) is practically con- 
stant and does not impose any further structure on the 
momentum distribution of the dissociation state. 



D. Connection to spectroscopy 

We now give an estimate of the matrix element 
(Prelj 0, O|bgW|0 res ) in terms of spectroscopically avail- 
able quantities such as the width of the Feshbach res- 
onance and the background channel scattering length. 
A natural basis of energy eigenstates is provided by the 
scattering states |0 P + '}, where [p 2 cl /2[i + Vb g (r-)]|0p + ^) = 

p 2 /(2/i)|</> p + ^) [22]. It describes the scattering in the rel- 
ative motion in absence of the confining laser poten- 
tials, with an incoming plane wave with momentum p 
as boundary condition. 

In order to relate (p re i, 0, 0| W\(p Tes ) to the spectroscop- 
ically available matrix element (0 p + ^|W|</> rcs ), we note 
that the spatial extension of the longitudinal trap is much 



larger than the support of the channel coupling W\(f> res ). 
This scale separation permits one to approximate the ef- 
fect of the longitudinal trap by a mere energy shift, 

(Prci,O,O| bg VF|0 ros ) pa (p rc i+PT,0,0| w/oT VF|</) res ), (23) 

where |p r ei, 0, 0) w / o x denotes the background channel en- 
ergy eigenstates of the relative motion without the trap 
potential (without loss of generality p lc \ > 0). The 
momentum shift pt is related to the trap depth by 
Uo,t = Inserting the identity in terms of the 

basis |</ip + '), we can rewrite the matrix element as 

res / 

d 3 P (Proi +PT,O,O| w/oT p + ))(0 p +)|IF|0 rcs ). (24) 

The matrix element (p rc i + pt, 0, 0| w / O T^p + ' > ) can be 
evaluated for a vanishing background channel potential, 
Vb g (r) = 0, since the scattering due to Vbg should not 
result in a substantial modification of the overlap. We 
thus have 



(Prol+PT,0, O| w / o T0p +) ) 

~ (Prel+PT\Pz){Pofl\Px)\Pv) 
1 

= 5(Prel +PT —Pz)- 



(25) 



-(pI+pD/2/j.Hwg 



Insertion into (24]) yields (p re i + pr, 0, 0| w/oT W\4> res ) « 

v 7 4irnhojG (</>q + ' ) | W | (pros) , where we use that p re i +px still 
lies in the ultracold regime such that W|<^ res ) cannot be 

\W\cf> ies ). This 



resolved and hence (</> P + ' ) | W\4>- 



is justified given that p Te \ is mainly determined by C. 
With (23) we thus obtain the desired connection to the 

^\W\ct> Ies ), 



spectroscopically available quantity 

CPrci,0,0|b g W^ ros ) ps V4vr^ G (^ +) |VF|</>r. 
According to Feshbach scattering theory [25] 

3 dbg /^res ^-^res 



\(4 +) \w\^ 



m(27rh) 



(26) 



(27) 



with ab g the background channel scattering length, 
AB rcs the resonance width, and /i ros the difference be- 
tween the magnetic moments of the Feshbach resonance 
state and a pair of asymptotically separated noninteract- 
ing atoms. 



E. Asymptotic relative motion 

Finally, let us approximate the basis states |p re i> 0, 0)b g 
in the asymptotic regime z vc \ — > oo, where the scattering 
states differ from the longitudinally free energy eigen- 
states only in a scattering phase (p S c{Prci), 



(aj re l|prel,0,0)bg ~ e 

z re \—>-oo 



^o(p r ol) (zrcl | 33rcl)(prcl |^rel ) _ 



(28) 
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The phase f S c(Prei) has two contributions, stemming 
from Vbg and Vt- Since we are in the ultracold regime of 
Vb g , its contribution is a linear shift given by the back- 
ground channel scattering length a^g. The contribution 
from Vt, on the other hand, can be linearized due to 
our requirement that the energy be sharply peaked, be- 
cause the width of the energy distribution is then small 
compared to the characteristic energy scale of the trap 
potential. The latter is determined by the trap depth 
Uq,t and is thus on the same order of magnitude as the 
kinetic energy after dissociation. This situation is simi- 
lar to the scattering of a narrow wave packet with spread 
AE at a broad resonance with width T, AE <§C T [3"4"j . 
Confinement-induced resonances [35] due to the guide 
potential are negligible provided the ground state size 
a± = (h/mus) 1 / 2 greatly exceeds the background chan- 
nel scattering length, a± ^> <2b g . 

As described by a linear scattering phase, the poten- 
tials thus merely effect an overall spatial displacement 
of the generated dissociation state. Physically, this shift 
stems from the faster propagation of the particles in the 
trap region. Since we are mainly interested in the struc- 
ture of the generated dissociation state, we can safely 
neglect this displacement and approximate 



(£Crel|Prel,0,0) 



bg 



(z IC l\p rC \)(Prcl\<P0 C fi 



(29) 



F. Canonical dissociation state 

Putting all together, the asymptotic dissociation state 
reads as 

|£Erel|— >00 III 
t—t-OO 

/CO poo 
dpcm / dp r d 
-OO J — OO 

x e -i(^bg+pL/2M+p r 2 ol /2 Al )t/?i 

x £ ^ + p§j2M+p^^ 

X (Pcm|V ; T)(^cmbcm)(Pcm|^S™) 
X (z rc l|prcl)(Prcl|^o C o)- (30) 

Normalizing the spectrum with 

/OO />oo 
dpcm / dp rcl \(p cm \ip T )\ 2 
-OO J —OO 

x \C (U hg /h + p 2 c j2Mh + p 2 j2[ih) | 2 (31) 
and introducing the abbreviation 

C bg = ^V47T M ?^G(4 +) |^l0re S }||(5||, (32) 

we can express the dissociation state in canonical form, 
I'M*)), « C hs uW\* z )\^ M%). (33) 

t— >oo 



Here U { °} = exp[-i(p 2 cm /2M + p 2 cl /2 t i + U hg )t/h] is the 
free longitudinal time evolution operator, and the longi- 
tudinal state component |?/> 2 ) is defined by the momen- 
tum representation 



(PcmKPrell^z 



C{U h Jh + plJ2Mh + p 2 j2^h) 



\C\\ 



X (PcmlV'l 



(34) 



The dissociation probability is given by |Cb g | 2 , which can 
be expressed in terms of the above mentioned spectro- 



scopic quantities using (27 1 



l^bgl = ^ \\c\\ ■ 



(35) 



It is mainly controlled by the applied magnetic field pulse 
which determines He'll 2 . We now focus on choices of the 
resonance width A_B ros and the magnetic field pulse such 
that the dissociation probability is on the order of a few 
percent. 

Taking Eqs. (33), (34 1 and (pBl together, we have 



found the desired expression for the asymptotic disso- 
ciation state. It is characterized by the trap geometry 
and, most importantly, by the Fourier transform of the 
closed-channel amplitude C(oj), which is in turn deter- 
mined by the applied magnetic field pulse sequence. In 
order to answer what kind of states can be generated, 
we thus have to determine the dynamics of the closed 
channel amplitude C(t). 



V. CLOSED-CHANNEL AMPLITUDE 
DYNAMICS 

In the previous section we found that the momentum 
representation ( 34 ) of the asymptotic dissociation state 



is mainly determined by the Fourier transform C(u>) of 
the closed-channel amplitude. In order to determine the 
generated state for a given magnetic field pulse sequence, 
we thus have to determine the dynamics of the closed- 
channel amplitude C(t) as determined by equation (16 1. 



A. Separation of decay and driving dynamics 



Let us rewrite the integral equation ( 16 ) as 
[ihd t ~ E ICS (B(t)) - U cl ] C(t) = f dt'f(t - t')C(t'), 

J —OO 

(36) 

with the kernel 



f(r) = -(^T|(^|(0rcs|W bg (r)lV-^T)|^!o>l^c S ) 



^E 1 



-iEr/h 



l(£|bgW^TMPI0rc S )| 



(37) 



where we used the decomposition of the background 
channel time evolution operator (17). The right-hand 
side of Eq. ( 36 ) describes the effect of the coupling be- 



tween the channels on the closed-channel amplitude. It 
is quadratic in the interchannel coupling W and hence, 
for a weak coupling and sufficiently short dissociation 
windows, its effect is expected to yield a small correc- 
tion to the unperturbed dynamics given by the left-hand 
side. Physically, we expect it to describe the decay of 
the closed-channel amplitude due to the escaping wave 
packet in the background channel. The kernel (37) may 



be viewed as the time-dependent overlap between the 
"initial state," V7|V'T)|Voo)l < ^res), and its evolved version 
Ubg{T)W\ipT)\ipQ a ^)\(j) ICS ) , which vanishes at large times 
due to the unbounded propagation in the background 
channel. It does not depend on the external magnetic 
field B{t). 



The kernel ( 37 ) is expected to drop off on a microscopic 
("memory") time scale t m , which can be roughly esti- 
mated from the spatial width Ax of the closed-channel 
bound state |^ rcs ). Denoting the corresponding momen- 
tum uncertainty by Ap, one obtains the drop-off time 
scale t m = mAx 2 /h from Apt m /m — Ax and the uncer- 
tainty relation. The spatial width of the closed-channel 
bound state \(j) les ) is on the order of the closed-channel 
scattering length, with typical values on the order of 
Aa; sa 100 clq. Taking the mass of lithium atoms ( 6 Li) one 
thus gets the estimate t m f=s 10 ns, which should be com- 
pared to the inverse decay rate of the resonance, which is 
much greater in our case. Given this shortness of t mi one 
might consider taking the limit t m — ¥ 0, which is equiva- 
lent to setting /(t) oc 6(t), but it will become clear later 
in this article that this approximation is too crude and 
cannot even qualitatively account for the correct decay 
behavior. 

In order to separate the anticipated decay from the 
unitary dynamics due to the left-hand side, we switch 
over to a "comoving frame" defined by 



C(t) = C (t)D(t), 



(38) 



where the uncoupled closed-channel amplitude Co(t) fol- 
lows by definition from 

[iHdt - E Ies (B(t)) - U cl ] C (t) = 0, (39) 

which implies 

C (t) = Co(*o) exp f-i J dt'[E Ies (B(t')) + U cl ]/h 

'(40) 



Applying the ansatz ( 38 1 on ( 36 ) and using ( 39 ) and ( 40 1 , 
one finds that the evolution of the coupling dynamics is 
governed by 

ihd t D{t) = I dt'D(t')f(t - t') (41) 

J — CO 

x exp (if dt"[E ies {B(t")) + U cl }/h^j . 



Since D(t) is driven only by the coupling between the two 
channels, we expect it to vary slowly for sufficiently small 
interchannel coupling W, such that it can be considered 
constant to good approximation on the time scale t m of 
nonvanishing kernel f(t — t'). This allows us to pull D(i) 
out of the integral, leading to 

d t D(t) = a(t)D(t), 

with the (in general complex) coupling coefficient 

ih 



a(t)=^ I dt'f(t-t') 



exp UJ dt"[E Ies (B(t")) + U cl }/h^j . (42) 



By writing 



a(t) = -T(t)/2 - iAE(t)/h, 



(43) 



we make explicit that the real and imaginary parts of a(t) 
describe the decay rate T(t) and an energy shift AE(t), 
respectively, as induced by the coupling between the two 
channels. 



B. Decay dynamics 



One can evaluate the coupling coefficient (42) further 



in the case of sufficiently smooth steering of the magnetic 
field, such that the resonance energy varies slowly on the 
scale of the drop-off time t m , 



d h 
-E Ies (t)t m <- : 



(44) 



In the vicinity of the resonance B Yes one can linearize, 
E ies {t) = /x»(B(t) - B Ies ), (45) 

leading to 



d „, . h 

U.0 L m /-*TCS 



(46) 



This assumption allows one to approximate the integral 
in the exponent of ( 42 ) as 



df" [£ ros (B(t"))+U c i]/h& [E ICS {B(t)) + U cl ]{t-t')/h, 

(47) 

such that we can rewrite the coupling coefficient ( 42 ) as 

1 r°° 

a (t) = ± dt'&(t-t')f(t-t') e' ; [^«(B(t))+c/d] (t-t')/» 

This can be read as the Fourier transform of the product 
of /(t) and the Heaviside step function 0(t), implying 



i 7(E Tes (B(t)) + U cl 



2h' 
1 

2^h 



h 



(49) 



3C:- 
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where V denotes the Cauchy principal value. Making use 
of (37), the Fourier transform f(u>) = dte lLUt f(t) of 
thclicrncl reads 

E 

(50) 

We thus find the decay rate according to ( 43 ) to be given 
by 

2tt 



and 



r(t) =y E 5 ( E re S (B(t)) + U cl - E) 
E 

x \{E\ hs W\^)\^l)\^ cs )\\ 
and the coupling-induced energy shift by 



(51) 



AE(t)=J2v 



1 



E ies (B(t)) + U c 



X |(£| bg W^TM n o}l0rc 



(52) 



Equation (51) shows that a nonvanishing decay rate is 
obtained only when E rcs (t) + U c \ matches a background 
channel energy eigenvalue. In particular, the gap in 
the spectrum between the dissociation threshold and the 
highest excited bound state explains why decay occurs 
only when the resonance energy lingers above the contin- 
uum threshold. (We always stay off-tuned from bound 
states of Vb g .) This also explains why the naive approx- 
imation for the kernel, /(r) ~ fo5{r), is not applicable; 
since the Fourier transform of S(t) is constant, it can- 
not distinguish energies above and below the continuum 
threshold and thus predicts an unphysical decay below 
the threshold. 



Note that (51), which coincides for constant magnetic 



field with the decay rate of the corresponding Feshbach 
scattering resonance, may also be viewed as a gener- 
alized version of Fermi's Golden rule, where the decay 
rate is determined by the instantaneous resonance en- 
ergy E TCS (B(t)), which in turn is externally controlled via 
the magnetic field B(t). This coincidence is not acciden- 
tal, of course, since the limit of a slowly varying magnetic 
field ( 46 ) admits both interpretations. The condition ( 46 ) 



quantifies the applicability of this approximation. 

Let us now use the results of the preceeding section by 
specifying the energy eigenbasis according to (20). As- 
suming again that the resonance state energy E Tes (B(t)) 
only sweeps over energies in the vicinity of the back- 
ground channel continuum threshold, remaining in the 
ultracold energy regime, below the first excited transver- 
sal state and off-tuned from the highest bound state of 
Vbg, we can write 

AE(t) =|( Prcl = 0,n r G el = 0,m r G e] = 0\ hg W\^- cs )\ 2 



cm 

oo J — oo 



dPrel |(Pcm|^T}| 2 



(53) 



x V 



E r , 



,(B(t)) + U cl -U hg -^ 



2M 



2p 



m ~\(Prel = 0,n r G el = 0,m r G el = 0\ hg W\^ cs }\ 2 



2tt 

/oo poo 
dpcm / dp rc l |(PcmV"r)| 2 (54) 
-OO J —OO 

x S(E ICS (B(t)) + U d - U hg ~ vlJlM - Pl 2 cl /2/i). 

We have omitted transitions into bound states of the lon- 
gitudinal trap; they are negligible given the pulse sweeps 
sufficiently fast over the corresponding energies. More- 
over, one can arrange the dissociation pulse B(t) such 
that the offset of the resonance state energy from the 
background channel continuum threshold greatly exceeds 
the trap state momentum uncertainty ct Pj t for most of the 
time, E res (B(t)) + U c \ — Uh g > o% T /2M. In that case 
the integrals are dominated by p IC \ 3> j) cm , and we can 
approximate 



KPcmlV'T)! 2 « S(Pc 



(55) 



This yields a vanishing energy shift, AE(t) — 0, since the 
principal value integration cancels. For the decay rate, 
on the other hand, we find 



r(t) 



2WGabgA*rcsAi? r , 



(56) 



2// 



V E Tes (B(t)) - 2C/ ,t + Swt/2 - hu> G 
xQ(E ICS (B(t)) - 2U , T + hw T /2 - Hu G ), 

where we substituted | (p rc i = 0, n^ 1 — 0, m G ' = 0| 
W\& cs )\ 2 
and U c i = 

we find that the decay rate ( 56 ) is nonvanishing only for 



= w G abg/XresAB res /7r, U hg = -2U ,g + 2frw G 
-2J7q,t + hwT/2 — 2C/ ,g + fioJG- As expected, 



magnetic field values that lift the resonance state above 
the (longitudinal) background channel continuum thresh- 
old. The offset — 2[/o,t+^ w t/2— Hujg in the step function 
gives the energy to be provided in addition to the free dis- 
sociation threshold; there, the trap depth to be overcome 
by both atoms is reduced by the closed-channel center 
of mass ground state energy Huj^/2, while the transver- 
sal relative motion, tightly bound in the closed channel, 
must make the transition to the transversal ground state 
of the guide. The square root pole stems from the one- 
dimensional state density and does not lead to apprecia- 
ble effects as long as the pulse sweeps sufficiently fast 
over it. 

In summary we find, under appropriate conditions on 
the dissociation pulse, that the decay dynamics of the 
closed-channel amplitude is described by 



d t D(t) 



r(t) 



D(t), 



(57) 



with the decay rate T(t) given by ( |56| ). By noting the 
formal solution D(t) = D(t )exp (- J" t * dtT(t')/2), the 
overall closed-channel amplitude C(t) then follows from 
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C(t) — Co(t)D(t), with the uncoupled closed channel 



amplitude Cb(t) given by (40) 



For the asymptotic dissociation state as described in 
Section |TV] we ultimately need to know the Fourier trans- 
form of the closed-channel amplitude, which is given by 
the convolution of the Fourier transforms of Co(t) and 
D(t), 



C{u>) = 



1 

2^ 



dujCo(tu)D(uj — oS). 



(58) 



In the limiting case of strong interchannel coupling or 
long-lasting, slowly varying dissociation pulses, the ki- 
netic energy distribution of the dissociated atoms, de- 
noted by n(E), is determined by the decay dynamics 
D (t) alone. In this quasi-stationary situation one may 
take the dissociation to occur monoenergetically, at the 
momentary resonance energy, by writing n (E) dE = 
r (t) \D (t)\ dt. For a monotonically increasing pulse en- 
ergy E TCS (t) the inverse t(E) exists (defining E Tes (E) = 
dtE les (t(E))) , and also the decay rate (56) can be viewed 



as a function of energy. Since D (t) decays exponentially 
it then follows that the energy distribution is given by 



n(E) 



dE 



exp 



Bo ^res (£') 



-dE' 



(59) 



This kind of quasi-stationary approach was used in [53] 
for the case of a linear field sweep, E rcs = const, and the 
above formula is consistent with their treatment when 
evaluated in the absence of confining lasers. 

On the other hand, in the case of a sudden magnetic 
field jump to a constant value Bq + AB above the thresh- 
old (which can be considered a square-shaped pulse in the 
limit of infinite pulse duration), the Fourier transform 
Co(w) gets sharply peaked at E ICS (Bo + AB), as will be 
shown later in this article. The function C(w) then re- 
duces to C(u) w D(uj - E ICS (B a + AB)/h), with D(u) 
a Lorentzian according to D(t) cx exp(— r(£' ros (_Bo + 
AB))t/2), which recovers the corresponding situation in 
[22]. 

However, in the following we are interested in the case 
where a magnetic field pulse dissociates on average only a 
single molecule out of the BEC. This implies that the de- 
cay of D(t) can be neglected compared to the dynamics of 
Co(t), such that the Fourier transform C(w) is essentially 
given by Co (w) . We will therefore focus on the uncoupled 
closed-channel amplitude Co(t), from now on, and in par- 
ticular on magnetic field pulses that result in a sharply 
peaked momentum distribution, as required for interfer- 
ometric purposes. The following section is devoted to 
magnetic field pulses that optimize Cq(uj) with that re- 
spect. Equation ( 58 1 shows that any non-negligible decay 



of the closed channel amplitude then merely results in an 
undesired smearing of the spectrum. 



VI. OPTIMAL MAGNETIC FIELD PULSE 

We proceed to characterize the magnetic field pulse 
shape that is optimal in terms of providing wave packets 
with well-defined momentum in the undepleted molecu- 
lar field limit. To this end, we can restrict the discussion 
to the case of a single dissociation pulse, since the gen- 
eralization to sequences of pulses is straightforward. In 
addition, as explained above, we are interested in situa- 
tions where the depletion of the BEC can be neglected, 
such that the closed-channel amplitude C(t) is well ap- 
proximated by the uncoupled component Co(i), whose 
dynamics is given by ( 39 ) and ( 40 ) . 



A. Dimensionless formulation 

We start by switching to a more convenient represen- 
tation. Since the resonance state energy E ICS (B(t)) is 
proportional to the magnetic field B(t) in the vicinity of 
the Feshbach resonance, see (45), it is sufficient to in- 
vestigate the dependence on E ICS (B(t)). We rewrite the 



integrand in the exponent of ( 40 ) as 



E ICS (B(t)) + U cl = E + AEP(t/T), (60) 

where P(t/T) describes a pulse with unit height; E de- 
notes the base value, and AE and T characterize the 
height and the width of the pulse, respectively. Intro- 
ducing the dimensionless energy e = AET/h and the 
phase function 



<f>(t/T) = f ' dtP(t)+0 o , 

Jt /T 

we can write the Fourier transform of Co (t) as 



(61) 



e -iE t /h _ roo 

- - C (uj + E /h)= / dte wt e~^ T \ (62) 

Given that the pulse function P(t) is positive and has 
a compact support, the phase function (f>(t) undergoes 
a monotonic ascent interpolating between two constant 
levels. If the pulse is also symmetric, P(—t) — P(t), the 
constant of integration (j)Q can be chosen such that the 
phase function (f>(t) is antisymmetric, 4>{—t) = —<f>(t). 



The constant prefactor on the left-hand side of ( 62 ) is 
irrelevant and will be neglected in the following. 



B. Gaussian magnetic field pulse 

We now ask for phase functions (f>(t) that yield a well- 
behaved spectrum Co(w). A natural starting point is to 
consider the spectrum resulting from a Gaussian-shaped 
pulse function P(t). In this case, as for most other pulses, 
the integral on the right-hand side of ( 62 ) cannot be cal- 
culated exactly. But asymptotic expansion techniques 
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can be applied if we take the dimensionless energy e 
to be large, £ > 1. This is a justified assumption in 
realistic scenarios. An estimate with ^ rcs w 10~ 2 /!b, 
AS = 100 mG (where B{t) = B + ABP(t/T)) and 
T = 100ms yields e — 10 3 . 

The Gaussian pulse leads to the error function, 

0(t) = erf(t), 

whose analyticity admits a uniform asymptotic expansion 
[3(3 ! which is capable of handling stationary points both 
lying on the real axis and in the complex plane. Dropping 
the prefactor on the left-hand side of (62 1, the result reads 



C (u + Eo/H) 



with the coefficients 



= 2m ^L M {e^W), (63) 



7 = 



i (| [erf (a) 
-(f[^«- 



uiT , 



<]) 

rfi(a)]) 



1/3 



1/3 



and 



where 



. (| [crf(a)-u;Ta/e]) 1/6 
' (uiTa/e) 1 / 2 

[wTQ/£-crfi(a)]) 1/6 



(wTa/e) 1 /2 



-ln(^ 



In 



if 



if 



if ^ < 2 

if ^ > 4= 



if «x < 

e — 
if !f > 

jT < _2_ 



> 



2 



Numerical analysis verifies excellent agreement with the 
exact result for large e. A plot of the corresponding spec- 
trum for e — 100 is given in Figure [2] It clearly does 
not exhibit a well-behaved peak structure, but rather 
shows contributions from all energies covered by the pulse 
sweep. The oscillating structure can be understood as 
due to the interference between the contributions from 
the ascending and the descending slope of the pulse. 



C. Square-shaped magnetic field pulse 

A clearer and more general insight into the relation 
between the pulse and the corresponding spectrum can 
be obtained by retreating to the stationary phase ap- 
proximation. This comes at the cost of losing the spec- 
trum in the tail region where the stationary points of 
the exponent in (62 1 leave the real axis. On the other 



hand, also non-analytic functions can now be treated. 
For the sake of simplicity, we consider symmetric pulses, 
P(—t) = P(t). In stationary phase approximation, we 



then get from ( 62 ) 



C (uj + E /h) 
T 



e\P'(t u )\ 



uTL, 



£<t>(t u 



for ojT/e < <^ 



(64) 




FIG. 2: Fourier transform Co{uj) of the uncoupled closed- 
channel amplitude, here evaluated for a Gaussian magnetic 
field pulse with the dimensionless energy e = 100. In the 
plot the numerically exact result is indis ting uishable from the 
asympt otic evaluation as given by Eq. ( 63 1 . As described in 
Section IV Cq{uj) essentially determines the momentum dis- 
tribution of the asymptotic dissociation state, see Eq. ( 34 1 . 



One observes that the spectrum does not exhibit a narrow 
peak-structure as one would require for atoms propagating 
with well-defined momenta. Rather, the spectrum exhibits 
contributions from all energies the pulse sweeps over. The 
oscillations stem from the interference between the contribu- 
tions of the ascending and the descending slope of the pulse. 



where t u denotes the positive (dimensionless) stationary 
point, defined implicitly by <p'(t u ) = ujT/e. Since t u is 
associated with the instant at which the pulse sweeps 
over the frequency w, it is manifest from (64 1 that the 



corresponding contribution to the spectrum is the more 
suppressed the faster the pulse sweeps over the corre- 
sponding energy, as one expects physically. 

Aiming at a spectrum that is well-peaked at a specific 
energy, we thus must strive for magnetic field pulses that 
sweep as fast as possible over the region of undesired 
energies and then rest at a plateau determined by the 
desired energy. The (idealized) optimal pulse with that 
respect is a square-shaped pulse, 



P(t) =e(t+l/2)0(t- 1/2), 



(65) 



for which ( 62 ) can be calculated directly, yielding 



C {uj + E /h) 
T 



—2tt cos 



smc 



£ 

2 



5{ujT) 



The pole at wT = can be traced back to the fact 



that the pulse function ( 65 ) is nonvanishing only on a 



finite time interval. For Eq < 0, which corresponds to an 
asymptotic magnetic field in the bound regime, this pole 
lies at a negative energy irrelevant for the shape of the 
dissociation state, as was shown in Section |IV| For our 
purposes it is thus sufficient to restrict the discussion to 
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FIG. 3: Fourier transform Co(w) of the uncoupled closed- 
channel amplitude component for a square-shaped magnetic 
field pulse with dimensionless energy e = 100. In contrast to 
the Gaussian pulse spectrum, it exhibits a pr onounced peak 
at uiT = e. Together with the connection (34 1 between C(uj) 



and the asymptotic dissociation state, this proves the capa- 
bility of the scheme to provide relatively well-defined atom 
momenta, even though there is no way to tailor arbitrary dis- 
sociation states. In particular, sequences of such pulses then 
can serve to generate highly nonclassical motionally entangled 
states. (Since Eq < the pole at ujT — does not affect the 
asymptotic behavior of the atoms.) 



positive frequencies with lj + Eq/H > 0. A plot of (66) 
for e = 100 is shown in Figure [3] As anticipated, the 
spectrum exhibits a pronounced peak at ljT = e, whose 
width is characterized by the pulse duration T. Accord- 
ing to (34), the corresponding dissociation state there- 
fore exhibits a narrow momentum distribution. This was 
used in [32] where the two-particle momentum distribu- 
tion generated by a square pulse has been worked out for 
a specific experimental scenario. 

This concludes the quest for the optimal magnetic field 
pulse shape. As indicated in the stationary phase calcu- 
lation ( [64] ) , by sweeping as fast as possible over undesired 
energies one singles out a spectrum peaked at the desired 
energy. This has been verified numerically for a large va- 
riety of pulse forms. In particular, it is now clear that 
smoothcning the edges of the square shape inevitably en- 
hances the undesired ripples in the region of the energies 
swept over, contradicting the naive expectation that they 
stem from the non-diffcrentiable ansatz ( 65 ) for the pulse. 

More generally, we found that the nature of Feshbach 
dissociation dynamics puts strong limitations on the pos- 
sible range of states that can be generated. It is in gen- 
eral not possible to find a pulse shape so as to generate 



a specified dissociation state. On the other hand, as was 
demonstrated in [25], the combination of a sequence of 
two dissociation pulses yields a promising perspective to 
generate highly nonclassical motionally entangled two- 
atom states with the capacity to violate a Bell inequality. 
VII. CONCLUSIONS 

We described the dissociation of Feshbach molecules 
as triggered by a time-varying external magnetic field in 
a realistic trap and guide setting. Employing a coupled 
channels formulation including the center of mass mo- 
tion, we derived an analytic expression for the asymp- 
totic dissociation state in the laser guide. It is essen- 
tially determined by the Fourier transform of the closed- 
channel amplitude C{t). The established relation be- 
tween applied dissociation pulse and resulting spectrum 
C{uj) puts strong limitations on the possibility to tailor a 
pulse shape in order to obtain a desired dissociation state. 
However, a square-shaped magnetic field pulse is shown 
to optimize the momentum distribution with respect to 
its width, as it is required for subsequent interferometric 
processing. 

The use of the Feshbach dissociation scheme as a 
source for entangled atom states was demonstrated in 
[25] . There, a sequence of two equal magnetic field pulses 
dissociates the molecules such that each atom is delocal- 
ized into a pair of macroscopically separated, consecutive 
wave packets. Their subsequent interferometric process- 
ing on each side then constitutes a Bell test whose success 
requires a sharp momentum distribution as provided by 
the discussed square-shaped pulses. 

An alternative protocol that can do without the in- 
terferometers employs a sequence of dissociation pulses 
of different heights. Given that the second pulse is the 
stronger one, the corresponding wave packets propagate 
faster as compared to the ones associated with the first 
pulse, such that time evolution alone effects their overlap. 
Position measurements in the overlap region on each side 
and their appropriate dichotomization may again reveal 
nonclassical correlations that entail a Bell violation. 

One may also discard the restriction to the dissociation 
of a single pair and consider the simultaneous dissociation 
of a multitude of molecules in a multipair Bell setting 
[57] . Statistical effects due to the indistinguishability of 
the particles then require a theoretical treatment beyond 
the single-molecule dissociation dynamics. For this, the 
above derived two-body transition amplitude (between 
molecular initial state and asymptotic dissociation state) 
could serve as an essential building block, for example 
with respect to a cumulant expansion [38] , 
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